function [H,indices]=dendrograma(k, g)
    %load iris.dat
    d = pdist(g,'euclidean'); % data son los datos
    %y = linkage(d,'average');  
    y = linkage(d, 'complete');
    set(0,'RecursionLimit',1500);
    [H,indices] = dendrogram(y,k*2); %, 'orientation', 'right'); 
